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Abstract. We investigate the interaction of a plane, linearly polarized, Alfvenic pulse with a one-dimensional, 
perpendicular to the magnetic field, plasma density inhomogeneity in the strongly nonlinear regime. Our numerical 
study of the full MHD equations shows that: (i) Plasma density inhomogeneity substantially enhances (by about 
a factor of 2) the generation of longitudinal compressive waves, (ii) Attained maximal values of the generated 
transverse compressive perturbations are insensitive to the strength of the plasma density inhomogeneity, plasma (3 
and the initial amplitude of the Alfven wave. Typically, they reach about 40% of the initial Alfven wave amplitude, 
(iii) Attained maximal values of the generated relative density perturbations are within 20-40% for 0.5 < j3 < 2.0. 
They depend upon plasma (3 strongly; and scale almost linearly with the initial Alfven wave amplitude. 
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1. Introduction 

The Alfven waves are usual candidates for energy trans- 
port from the lower layers of the solar a t mosp here to the 

~ 2000| ). However, 



• corona, e.g. (Goossens 1994; Roberts 



efficient deposition of the momentum and energy re- 
quire interaction of linearly incompressible Alfven waves 



with compressible magnetoacoustic waves, e.g. (Ofman & 
Davila 1 |1997|, |1998|; pfman et aT~| |2000|). Also, the com- 



pressible waves, in contrast to the Alfven waves, can trans- 
port energy and momentum across the magnetic field, 
spreading out the heated region. In addition, observational 
detecti on of Alfven wavoo in open otrueturoo of the corona 
can bel baocd upon mcaouromont of the comprcooiblo flue 



tuations, e.g. (Ofman ct al. 1997, 1998, 2000). These 



are generated in the lower corona by the Alfven waves 



through linear or nonlinear mechanisms, e.g. (Nakariakov 



et al. 2000). This method can be complimentary to the 



observation of coronal Alfven wav es through non-therm al 



broadening of emission lines, e.g. ( |Banerjee et al 



In the inner heliospheric solar wind, Alfven waves are 
observed in situ and represent the main component in 



MHD turbulence (Tu fc Marsch 1995; Tsurutani fc Ho 



1999). This suggests another interesting problem: why 



the incompressible turbulence dominates in the solar wind, 
and why compressible fluctuations are not observed, de- 
spite the theoretical possibility for these two kinds of MHD 
fluctuations to be coupled because of the medium inhomo- 
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geneity and nonlinearity. Thus, the study of coupling of 
compressible and incompressible fluctuations is important 
to the physics of the solar corona and the solar wind. There 
are several possible mechanisms for the coupling. 

The decay instability of Alfven waves is one of the 
possible examples of interaction between the MHD wave 
modes. This mechanism involves resonant three-wave 
interaction of Alfven and magnetoacoustic waves (e.g. 
[Sagdeev fc Galeev" ( 1969| ) and references therein). The ef- 
ficiency of this interaction is governed by the amplitudes 
of the interacting waves. However, this mechan ism works 
only for qu asi-periodic (perhaps wide-spectrum, Malara et 



al 



(200C)) waves, and is not efficient for the wave pulses 
that could be generated by some transient events on the 
Sun such as solar flares and coronal mass ejections, e.g. 



( [Roussev et al. 2001 ) 
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In contrast, the efficiency of non-resonant mechanisms 
of the compressible fluctuation excitation by the Alfven 
waves does not depend on the coherentness of Alfven per- 
turbation and consequently work even for single wave, 
wide-spectrum, pulses. 

Nonlinear excitation of magnetoacoustic perturbations 
by nonlinear elliptically polarized Alfven waves via the 
longitudinal gradients of the total pressure perturbations 
is one of the possible examples of the non-resonant MHD 
wave interaction. In this mechanism the generation of 
compressible perturbations results in the self-interaction 
and subsequent steepening of the Alfven wave front, which 
is described by the Cohcn-Kulsrud equatio n, e.g. (pohen 
& Kulsrud 1 119741; IVerwichte et al. 1 11999|; iNakariakov et 



2 



Tsiklauri, Nakariakov, & Arber: A Strongly Nonlinear Alfvenic Pulse 



al. 2000). In the following, we refer to this mechanism 
as "longitudinal", to distinguish it from the "transverse 
mechanism" , which generates compressible perturbations 
via the transverse gradients of the total pressure perturba- 
tions in nonlinear Alfven waves. More detailed discussion 



of these two mechanisms is presented in Nakariakov et al. 
( |1997| , [1998D ; [Botha et al. | ( |2000| ); [Tsiklauri ct al. | ( |2001| ) 



The transverse mechanism is dramatically modified in 
the case when MHD waves interact with transverse sta- 
tionary inhomogeneity of the plasma. If the Alfven speed 
is inhomogeneous across the magnetic field, initially plane, 
linearly polarized Alfven waves become oblique and sharp 
gradients in the direction across the field are secularly gen- 
erated. This constitutes a basis of the well known Alfven 
wave phase mixing phenomenon which has been exten- 
sively investigated in connection with the solar coronal 
heating with MHD waves ( |Heyvaerts fc Priest | |1983| ). 



Various aspects of this phenomenon have been intensively 



st udied using full MHD numerical simulations ( Malara et 
al. ||1906|; |Ofman & Davila ||l995|, |1997|; |Pocdts et al.||l997|; 



Grappin et al. 


200C|; Dc Moortcl ct al. 


2000). As demon- 


strated by 


Nakariakov et al. 


(|1997 


, 1998|); Botha et al. 



(|2000|) ; [Tsiklauri et al. | (J2001J) , phase mixing of Alfven 
waves in the compressible plasma, in the weakly nonlin- 
ear regime leads to the generation of fast magnetoacoustic 
waves, and various regimes of this process, relevant to solar 
coronal and heliospheric applications have been studied. 
In particular, it has been found that the inhomogeneity of 
the plasma across the magnetic field, associated with var- 
ious types of structuring (e.g. plumes in the coronal holes, 
boundaries between the slow and the fast solar winds, flow 
tubes, etc.), plays the crucial role in the interaction of 
compressible and incompressible weakly-nonlinear MHD 
modes. 

Problems connected with the interpretation of MHD 
fluctuations observed in the solar wind and the propaga- 
tion of intensive MHD wave pulses in the solar corona 
(e.g., flare and CME-associated waves) require detailed 
study of the interaction of the MHD pulses with plasma 
inhomogeneities in the strongly nonlinear regime. Our in- 
terest to this regime is motivated by the necessity to ac- 
count for effects of higher order nonlinearities and the 
back-reaction of nonlinearly generated compressible per- 
turbations on the source Alfven pulse. Also, as the effi- 
ciency of the "transverse" mechanism is connected to the 
longitudinal wave number (together with the transverse 
wave number and the amplitude), and as the "longitudi- 
nal" mechanism increases, through wave steepening, the 
longitudinal wave numbers, we expect that the simulta- 
neous action of these two mechanisms can enhance the 
efficiency of the transverse mechanism. In this work, we 
study, by direct numerical simulations, the generation of 
compressible fluctuations by a strongly nonlinear Alfvenic 
pulse interacting with a transverse plasma inhomogeneity. 

The paper is organized as follows: in Section 2 we de- 
scribe our model and the numerical method applied, in 
Section 3 the results of the simulations are discussed sep- 
arately in the high and low (3 cases (subsections 3.1 and 



3.2 respectively), sub-section 3.3 deals with the investiga- 
tion of parametric space, and finally, the conclusions are 
presented in Section 4. 

2. The model 

The model studied here is similar to one discussed in 
Tsiklauri et al] ( |200 1| ): we use the equations of ideal MHD 



dV 
} ~dt 



p(V ■ V) V = —Vp — -— -B x curLB, (1) 



47T 



dB 

~dt 



curl(V x B), 



7pV • V = 0, 



^+div( P V) = 0, 



(2) 
(3) 
(4) 



where B is the magnetic field, V is plasma velocity, p is 
plasma mass density, and p is plasma thermal pressure for 
which the adiabatic variation law is assumed. 

We solve equations (l)-(4) in Cartesian coordinates 
(x, y, z) and under the assumption that there is no varia- 
tion of the physical values in the y-direction, i.e. (d/dy = 



0, 2.5D approximation) with the use of Lare2d (Arber et 
|al. |2001 ). Lare2d is a numerical code which operates by 
taking a Lagrangian predictor-corrector time step and af- 
ter each Lagrangian step all variables are conservatively 
re-mapped back onto the original Eulerian grid using Van 
Leer gradient limitcrs. This code was also used to produce 

~ (|200q);|Tsiklauri et al 



the results in Rcfs. 



(2001). As in Botha et al 



Botha et al 

pjoo[ ) 



Tsiklauri et al. (2001), 



the equilibrium state is taken to be a uniform magnetic 
field Bq in the z-direction and an inhomogeneous plasma 
of density po(x), 



Po(x) = p» (3 — 2 tanh(Aa;)) 



(5) 



Here, A is a free parameter which controls the steep- 
ness of the density profile gradient. The latter is lo- 
calized around x = 0. The temperature profile Tq(x) 
is set up to allow for the total pressure to be con- 
stant everywhere. In our normalization, which is the 
same as that of |Botha ct al. | ( |2000[ ); [Tsiklauri ct] 
FH (pOOlp, B = BqB, r = 
Bo/y/inp (x) = [Bo/^AWp-;] / ^3 



C A (x) 



C" 



[Bo/^/WpZ 
3 - 2 tanh(Ax) = C* A C A {x), t 



- 2 tanh(Ax) 
(a*/C* A )i, V 



C A V, the dimensionless local Alfven speed is Ca(x) = 
l/y/Z-2 tanh(Ax). a* and p* are the units of length and 
density respectively. In what follows we omit bars on all 
dimensionless physical quantities. 



3. Numerical Results 

We set up the code in such a way that initially longitudi- 
nal (V z ) and transverse (V X ,B X ,B Z ) perturbations and the 
density perturbation, p, (perturbed by both longitudinal 
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and transverse modes) are absent and the initial ampli- 
tude of the Alfven pulse is strongly non-linear, i.e. typi- 
cally A=0.5. At t = 0, the Alfven perturbation is a plane 
(with respect to x-coordinate) pulse, which has a Gaussian 
structure in the z-coordinate, 



By(z) = Aexp ( -— ) , V y (x, z) = -C A (x) B y (z). (6) 



Here, S is a free parameter which controls the width of the 
initial pulse. In all our numerical runs S = 0.1. 



As it is discussed in Introduction (cf. Nakariakov et 



al. (It 97) for details) in the considered geometry, Alfven 



waves represented by V y and B y , are incompressible up to 
the cubic nonlinearity. The longitudinal (represented by 
V z ) and the transverse (represented by V X ,B X ,B Z ) com- 
pressible perturbations (of course, both perturbing p) are 
generated nonlinearly by the longitudinal and transverse 
gradients of the total pressure. The efficiency of the gener- 
ation depends upon the amplitude of the Alfven waves, in 
other words, the incompressible and compressible pertur- 
bations are linearly decoupled. The latter guarantees that 
with the choice of our initial conditions the compressible 
perturbations are indeed initially absent from the sys- 
tem, which a priori is not clear if the linear coupling is 
present (cf. [Tsiklauri fc Nakariakov ( [2002 )). In addition, 
the transverse compressible perturbations are generated 
by plane Alfven waves only in the presence of the trans- 
verse profile of the local Alfven speed. 

The simulation box size is set by the limits —15.0 < 
x < 15.0 and —15.0 < z < 15.0. The pulse starts to move 
from point z = —12.5 towards the positive z's. Both the 
scale of density inhomogeneity and width of the Alfvenic 
pulse are much smaller than the size of the calculation 
domain. 

We have performed calculations on various resolutions 
in an attempt to achieve convergence of the results. The 
graphical results presented here are for the spatial resolu- 
tion 2500 x 2500, which refers to number of grid points in z 
and x directions respectively. We have used a non-uniform 
grid in our simulations, namely, in x direction 75% of the 
grid points where concentrated between —5.0 < x < 5.0 
where the spatial inhomogeneity has strongest gradients. 
We have also performed calculation on the spatial resolu- 
tion 3500 x 3500 and found that the maximal generation 
levels (maximum of an absolute value over the whole sim- 
ulation domain) for all physical quantities as a function of 
time are the same as in the case of 2500 x 2500 resolution. 
Thus, the results presented here are, indeed, converged. 

Our main numerical runs are split in two parts, treat- 
ing separate cases when plasma (3, which is the ratio of 
speed of sound to the Alfven velocity squared, is less 
and greater than unity. This split is motivated by weakly 
nonlinear results, showing that the particular scenario of 
the interaction of MHD waves is determined by the ratio 



(C s /Ca) being less or greater then unity, e.g. (Cohen & 
Kulsru|^ |i974l ). 



Fig. 1. Top panel: snapshot of V x /Ca(x) at t — 15.0. 
Bottom panel: contour-plot of V x /Ca{x) at the same in- 
stance. Here, initial amplitude, A = 0.5, density inhomo- 
geneity steepness, A = 0.75, plasma (3 = 2. 

Fig. 2. Top panel: snapshot of V z /Ca(x) at t = 15.0. 
Bottom panel: contour-plot of V z /Ca(x) at the same in- 
stance. Here, the parameters are the same is in Fig. 1. 

Fig. 3. Top panel: snapshot of (p — po)/po at i = 15.0. 
Bottom panel: contour-plot of (p — po)/po at the same 
instance. Here, the parameters are the same is in Fig. 1. 



3.1. Case when [3 > 1 

In this subsection we present solution of the Eqs. (l)-(4) 
with the above described equilibrium and the initial condi- 
tions for the case when plasma-/3 is 2.0. Here, A-parameter 
was fixed at 0.75. 

Figures 1-3 show snapshots of the initially absent 
transverse (V x ) and longitudinal (V z ) components of the 
velocity (representing the "transverse" and "longitudinal" 
compressible waves, respectively) and density perturba- 
tion, at time £ = 15. 

The transverse compressible perturbations (V x ) are 
generated by the transverse gradients of the total pres- 
sure perturbations. In turn, these perturbations are gen- 
erated by the Alfven wave phase mixing, and are located 
near x — where the density gradients are the largest. 
Then they propagate across the field (see Fig. 1). The 
longitudinal compressive perturbations, V z are generated 
even in the absence of the plasma density inhomogeneity 
(Fig. 2). In the contour plot (Fig. 2, bottom panel) it can 
be seen that there are two wave fronts: one moving at a 
local Alfven speed, and another (compressive one) that 
moves faster than the local Alfven speed because (3 = 2. 
Fig. 3 shows the relative density perturbation associated 
with the compressible waves that have been generated. 

In Fig. 4, the top panel presents spatial variation 
(across the x-coordinate) and the dynamics of the trans- 
verse and longitudinal compressible waves, and relative 
density perturbations produced by them. In what follows 
max and min refer to the maximum and minimum over 
the whole simulation box (i.e. space) of a dimensionless 
physical quantity at a given time instance, respectively. 
These are fairly good and simple (scalar) quantities de- 
scribing the generation and/or decay of a 2.5D physical 
quantity, dynamics of which is otherwise not so straight- 
forward to comprehend. There are five interesting observa- 
tions: (i) the decay of the max{\V y \ / C a{x)) (on expense of 
which the transverse (V x ) and longitudinal (V z ) compres- 
sive waves and associated density perturbation are gen- 
erated) occurs in the middle - where the plasma inhomo- 
geneity is the strongest. This demonstrates the importance 
of the inhomogeneity. (ii) the right wing (x > 0) decays 
faster than the left (x < 0) one. For the value of plasma (3 
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used, it is expected that shock dissipation would be greater 
where the local Alfven velocity is greater (in this case the 
right wing). Note that when referring to shock dissipa- 
tion we mean artificial dissipation that guarantees proper 
shock capturing, while at all times we remain in the frame- 
work of an ideal MHD (no bulk dissipation included). 
Thus, the shock viscosity ensures that we recover the weak 
solution to ideal MHD. (iii) the relative density perturba- 
tion and longitudinal compressive wave max(\V z \/CA(x)) 
are generated in the homogeneous parts of the domain 
(x < —5, x > 5) too. However, in the middle of the do- 
main, where the inhomogcncity is strong we see that in- 
homogeneity of plasma density enhances the generation of 
these quantities by about a factor of 2. (iv) transverse com- 
pressive wave max(\V x \/CA(x)), which is not generated in 
the absence of density inhomogcncity, is generated in the 
middle of the domain, (v) max(\V y \/CA(x)), on expense of 
which the transverse (max(\V x \ / Ca(x))) and longitudinal 
{max{ | V z \/Ca{x))) compressive waves are generated, has 
two dips in those places where the other physical quan- 
tities have two bumps, which clearly demonstrates the 
correct energy balance as well as the importance of the 
inhomogeneity (where actually the bumps occur). 

In Fig. 4, the bottom panel presents temporal varia- 
tion of transverse and longitudinal compressible wave am- 
plitudes. In the presence of the inhomogeneity (A = 0.75), 
the Alfven perturbation (min(\V y (x, z,t)\/CA(x))) decays 
much faster than in the homogeneous case (A = 0). In 
the inhomogeneous case, the energy initially stored in the 
Alfven wave in addition to shock dissipation goes into the 
generation of the compressive waves. The maximum am- 
plitudes of longitudinal max(\V z (x, z, t)\) and transverse 
max{\V x {x, z,t)\) compressible waves attain a substantial 
fraction of the initial Alfven wave amplitude. In the inho- 
mogeneous case (A = 0.75), the longitudinal wave attains 
about twice the maximal value than in the homogeneous 
plasma case. Obviously, when A = there is no generation 
of the transverse compressive wave and max(\V x (x, z, t)\) 
is identically zero for all times. 

3.2. Case when (3 < 1 

In this subsection we present solution of the Eqs. (l)-(4) 
with the above described equilibrium and the initial condi- 
tions for the case when plasma-/? is 0.5. As in the previous 
subsection, A-parameter was fixed at 0.75. 

Again, the initially absent transverse (V x ) and longi- 
tudinal (V z ) compressive waves and density perturbations 
are efficiently generated (Figures 5-7) and reach a sub- 
stantial fraction of the initial Alfven wave amplitude. 

In the contour plots in Figs. 6 and 7 it can be seen that 
there are two wave fronts: one moving at a local Alfven 
speed, and another (compressive one) that moves slower 
than the local Alfven speed because (3 = 0.5 (compare 
with Figs. 2 and 3). Therefore, since the velocity difference 
between the first and second parts of the solution is not 



as great as in the /3 = 2 case, these two parts seem be 
blended into each other. 

Figure 8 presents the evolution of transverse and longi- 
tudinal compressive wave amplitudes in time. This graph 
is quite similar to Fig. 4. The noteworthy difference is as 
follows: for the value of plasma (3 used here {j3 = 0.5), the 
non-linear term in the scalar Cohen-Kulsrud equation is 
larger than in the former case of = 2 (Fig. 4) . Therefore 
in Fig. 8, top panel, we observe that the shock dissipation 
of the Alfven wave is greater (the thin solid line goes fur- 
ther down than in Fig. 4, also the top panel), and, in turn, 
the non-linear generation of the transverse and longitudi- 
nal compressive is further enhanced. 

3.3. Parametric study 

In this subsection we explore the parametric space of the 
problem. In particular, we investigate how the maximal 
value of the generated transverse compressive wave de- 
pends on the plasma density inhomogeneity parameter, A, 
plasma (3, and initial amplitude of the Alfven wave, A. 

In Fig. 9A we plot the dependence of the maximum 
of the absolute value of the transverse compressive per- 
turbation, max(\V x (x,z,t)\)/A, versus A for j3 — 2 (solid 
curve) and (3 = 0.5 (dashed curve). There are two note- 
worthy features in this graph, first, the maximal value of 
the generated transverse compressive wave depends on the 
plasma density inhomogeneity parameter rather weakly 
(once A > 0.4), and second, efficiency of the generation of 
V x is somewhat larger in the (3 = 0.5 case than in the case 
of f3 = 2. There are no data points between < A < 0.4 
because in order to accommodate the inhomogeneity in 
the simulation domain we would have to increase its size, 
which was not possible with available computational re- 
sources. Also, we made 5 runs of Lare2d code for the dif- 
ferent values of (3. The results are presented in Fig. 9B. We 
gather from the graph that the maximal value of the gen- 
erated transverse compressive wave depends on the plasma 
(3 rather weakly. 

Yet another valuable insight can be obtained by study- 
ing the dependence of the maximal value of generated 
transverse compressive wave on the initial amplitude of 
the Alfven wave as our problem is essentially non-linear. 
In Fig. 9C we plot results of numerical runs for differ- 
ent values of A, while A and (3 where fixed at 0.5 and 
1.25 respectively. We gather from Fig. 9C that quite un- 
expectedly the ratio max(\V x (x, z, t)\)/A is insensitive to 
the variation of the initial amplitude of the Alfven wave. 

We also investigated the parametric space with regard 
to the relative density perturbation. Namely, we investi- 
gate dependence of max(\(p — po) / po\) as function of A, (3 
and A. We gather from Fig. 9D that for (3 = 2 the gen- 
erated relative density perturbation saturates at about 20 
%, while for (3 = 0.5 the saturation level doubles. Fig. 9E 
illustrates the fact that the maximum generated relative 
density perturbation depends quite strongly on plasma (3. 
In Fig. 9F we show the dependence of max(\(p — po)/po\) 
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Fig. 9. (A) Dependence of max(\V x (x, z, t)\)/A versus A (density inhomogeneity steepness) for f3 = 2 (solid curve) and 
(3 = 0.5 (dashed curve). Here, initial amplitude A — 0.5. (B) The same, but versus (3 for A = 0.5 and A = 0.5. (C) The 
same, but versus A for A = 0.5 and (3 = 1.25. (D) Dependence of max(\(p — po)/po\) versus A for (3 — 2 (solid curve) 
and (3 = 0.5 (dashed curve) both for A = 0.5. (E) The same, but versus (3 for A = 0.5 and A = 0.5. (F) The same, but 
versus A for A = 0.5 and (3 = 1.25. 



versus the initial amplitude A. We observe that since the 
density perturbation is generated by the non-linear effects, 
max(\(p — po)/po\) indeed grows with the increase of A, 
and the dependence is almost linear. 



4. Conclusions 

This s tudy is an extension of the previous works to the 



case o j strongly nonlinear amplitudes. The main results 
can be summarized as follows: 

— Phase mixing of a strongly nonlinear Alfven pulse is 
accompanied by an enhanced generation of compress- 
ible waves. This is true irrespective of plasma (3 being 
less or greater than unity. 

— Plasma density inhomogeneity (while providing, as in 
the weakly nonlinear case, a source region for the non- 
linear generation of transverse compressive waves by a 
plane Alfven wave) substantially enhances (by about a 
factor of 2) the generation of longitudinal compressive 
waves. 

— Attained maximal values of the generated trans- 
verse compressive perturbations are insensitive to the 
strength of the plasma density inhomogeneity, plasma 
(3 and the initial amplitude of the Alfven wave. 
Typically, they reach about 40% of the initial Alfven 
wave amplitude. 



— Attained maximal values of the generated relative den- 
sity perturbations are within 20-40% for 0.5 < (3 < 2.0. 
They depend upon plasma (3 strongly; and scale (in- 
crease) almost linearly with the initial Alfven wave am- 
plitude. 

Commenting upon the maximal values attained by the 
transverse compressive waves in relation to the previous, 



weakly nonlinear, results (Botha et al. 2000; Tsiklaur: 



et al. 2001) we would like to state the following: In 



the weakly nonlinear case the transverse compressive 
waves do not grow to a substantial fraction of the initial 
Alfven wave amplitude due to the destructive wave in- 
terference effect (typically, max(\V x (x,z,t)\)/A ~ 10~ 3 ). 
However, in the strongly nonlinear case, studied here, al- 
though the destructive wave interference is still in action, 
rnax(\V x (x,z,t)\)/A is now significantly larger, about 0.4 
(cf. Figs. 9A, 9B, 9C). Therefore, we conclude that the 
strong nonlinearity substantially enhances the attained 
maximal values of the transverse compressive waves. 

The physical phenomenon studied here is an elemen- 
tary process responsible for non-resonant coupling of com- 
pressible and incompressible MHD modes. In particular, 
it may play a role in MHD turbulence of space and as- 
trophysical plasmas. The presence of a pressure-balanced 
inhomogeneity should significantly affect the saturated 
MHD turbulent state (cf. a similar claim but based upon 
different reasonings in Bhattacharjee et al. ( 1998| )). 
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Fig. 4. Top panel: spatial variation (across x-coordinate) 
and the dynamics of non-linear generation of the 
transverse and longitudinal compressive waves as well 
as relative density perturbation in time. Thick solid 
line corresponds to the initial value (at t = 0) 
of max(\V y \/CA(x)). Thin solid curve represents the 
same, but for t = 7.5. The relative density perturba- 
tion max(\(p — po)/po\) (dotted curve), the longitudinal 
compressive wave max(\V z \/Ca(x)) (dash-dotted curve), 
and the transverse compressive wave max(\V x \/Ca(x)) 
(dashed curve) are given for t = 7.5. Bottom panel: 
Evolution of max(\ V x (x, z, i)|), min(\V y (x, z, £)|/CU(aO), 
max(\V z (x, z, t)\) in time. The solid curve with open 
rectangles presents decay of the initial Alfven perturba- 
tion (min(\V y (x, z,t)\/CA(x))) due to shock dissipation 
in the case of absence of the plasma density inhomogene- 
ity (A = 0) for three time instances (t — 0,7.5,15.0). 
The dotted curve with open triangles depicts the same, 
but when A = 0.75. The dash-dotted curve represents 
max(\V z (x, z, i)|), while thick solid curve corresponds to 
max(\V x (x, z, t)\) both for the case of A = 0.75. The 
long dashed curve represents max(\V z (x, z, t)\) for the case 
when A = 0. Here, the parameters are the same is in Fig. 1. 
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Fig. 5. Top panel: snapshot of V x /Ca(x) at t = 15.0. 
Bottom panel: contour-plot of V x /Ca{x) at the same in- 
stance. Here, initial amplitude, A — 0.5, density inhomo- 
geneity steepness, A = 0.75, plasma [5 — 0.5. 



Fig. 6. Top panel: snapshot of V z /Ca{x) at t = 15.0. 
Bottom panel: contour-plot of V z /Ca{%) at the same in- 
stance. Here, the parameters are the same is in Fig. 5. 



Fig. 7. Top panel: snapshot of (p — po)/po at t = 15.0. 
Bottom panel: contour plot of (p — po)/po at the same 
instance. Here, the parameters are the same is in Fig. 5. 





Fig. 8. The same as in Fig. 4, but with (3 = 0.5. 



